On the mechanism for breaks in the cosmic ray spectrum 
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The proof of cosmic ray (CR) origin in supernova remnants (SNR) must hinge on full consistency of the CR 
acceleration theory with the observations; direct proof is impossible because of the orbit stochasticity of CR 
particles. Recent observations of a number of galactic SNR strongly support the SNR-CR connection in general 
and the Fermi mechanism of CR acceleration, in particular. However, many SNR expand into weakly ionized 

C^l ' dense gases, and so a significant revision of the mechanism is required to fit the data. We argue that strong 

' ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated 

particles by exactly one power. The spectral break is caused by a partial evanescence of Alfven waves that 
confine particles to the accelerator. The gamma-ray spectrum generated in collisions of the accelerated protons 
with the ambient gas is also calculated. Using the recent Fermi spacecraft observation of the SNR W44 as an 
example, we demonstrate that the parent proton spectrum is a classical test particle power law <x E~ 2 , steepening 

1—5 . to£- 3 atE br ^IGeV. 
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i 1 ■ I. INTRODUCTION 

<\ 

The discovery of cosmic rays (CR) dates back to the historic Victor Hess balloon ascent in 1912 |Q]]. CR origin is thus a 
century old problem. Only the latest direct observations of galactic supernova remnants (SNRs) |0-|9l narrowed the search to 
precisely these objects as the most probable sources of the CRs. One serious problem on the observational side was the lack of 
the SNR gamma-ray data below the energy range of the imaging atmospheric Cerenkov telescopes, or IACT. The Fermi-LAT 
(large area telescope) and Agile observatories are rapidly bridging this gap (roughly in the 0.1-30 GeV band e.g., fllill), virtually 
overlapping with the IACT energy band. There have been recent breakthrough observations of such SNR as W44, IC443, W28, 
RX J1713 and Cas A S [T3]. Overall, observations favor the diffusive shock acceleration (or DSA |fllHl3ll , a modern 
i— i ' version of the mechanism originally suggested by Fermi in 1949 O) as a means for the production of galactic CRs. However, 
there are questions, and even some challenges, that the recent observations pose to the theory. 

Of those, the most relevant to the proof of the SNR-CR connection is the form of the spectrum that the theory predicts for 
the particular SNR conditions. Full understanding of the spectra will allow one to disentangle the proton (i.e., the primary CR 
qq • component) emission from a contaminating (1-2% level) but radiatively more efficient, and accessible to the direct observations, 
£f~) ' electron CR component. The most recent challenge to the DSA was posed by the measurements of the rigidity (momentum 
t— I , to charge) spectra of different species (most notably proton and helium). They turned out to be different, contrary to the DSA 
' predictions for the ultra-relativistic rigidity range. 

Note that the latter problem arose from the indirect observations of the background CRs Il5l - [l8tl . as opposed to the above 
£N1 1 mentioned direct observations of the putative accelerators (SNR). Generally, it is impossible to trace CR back to their accel- 
erators because of the orbit scrambling. The proof of their origin in SNRs can only be achieved by proving the acceleration 
theory consistent with all accessible observations. It should be noted that 'direct' observations also provide only the secondary 
photon emission generated by accelerated particles, either electrons (through synchrotron, Bremsstrahlung and inverse Compton 
radiation), or protons (through their collisions with the ambient gas material). Therefore, such observations cannot be interpreted 
as an evidence of proton acceleration in SNR without a detailed understanding of the emission mechanism. Note that electron 
acceleration in SNRs to at least ^100 TeV is held proven "beyond a reasonable doubt" after the observations of the SNR 1006 
by ASCA and other X-ray instruments lfl9ll20ll . 

This paper deals with the modification of the DSA proton spectra in a partially ionized SNR environment and its signatures 
in gamma-emission from such remnants. The recent discovery of the proton/helium anomaly in the background CR spectra 
is discussed elsewhere ( 1I2II1 . see also 1I22I - I24I1 for other suggestions to explain this anomaly). Here we pursue an alternative, 
complementary approach to more common multi-band treatments, e.g., 02511 . where the fits are primarily focused on the overall 
agreement across the entire spectrum (from radio to gamma). By contrast, we concentrate on the gamma-ray band and fit an 
important signature of the spectrum which is the spectral break. We believe it conveys an important information about the physics 
of acceleration missed in the 'standard' DSA theory. The quality of our fit, with virtually no adjustable parameters, should testify 
for the underlying physical scenario behind the emission. The broad-band fits do not typically meet high-quality criteria, as they 
seek to fit several portions of the data simultaneously by adjusting, in some cases, a few free parameters. Nevertheless, they 
provide an excellent consistency check for each particular model. 

The recent Fermi-LAT observations of the SNRs W44 and IC443 10 13] indicate that the spectrum of the gamma ray producing 
protons is substantially steeper in its high energy part than the DSA predicts. A similar discrepancy has been found in the 
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high energy gamma ray spectra measured by e.g., the CANGAROO 0], HESS fif] and MAGIC l26ll atmospheric Cerenkov 
telescopes. The lack of understanding of the primary particle spectra triggered debates about the nature of the observed gamma- 
ray emission (hadronic vs leptonic), e.g., rf27ll . We argued lEH that when a SNR interacts with a dense molecular cloud complex, 
the conditions for particle confinement to the shock are different from those adopted in conventional DSA modeling. Since the 
propagation of resonant Alfven waves is inhibited by ion-neutral collisions, particles are not confined and so escape the emission 
volume. These phenomena should result in a spectral break in the parent proton and thus, in the gamma-ray spectrum. The 
spectral index at the break should change by exactly one power Aq — 1 due to an effective reduction of particle momentum 
space dimension by one, since particles are confined in coordinate space only when they are within a slab in momentum space 
oriented perpendicular to the local mean magnetic field. Note that the earlier HESS observations of the SNR RXJ 1713 were 
also consistent with such a break yfl. The most convincing evidence for the breaks of index one, however, provide the recent 
Fermi-LAT and Agile observations of W44 fiHIH (re-analyzed in (Q), the MAGIC observations of the SNR W51C l26l,|3lll 
as well as the FERMI observations j32tl of giant molecular clouds (GMC), where the Alfven wave evanescence should also 
result in a Aq ~ 1 steepening of the E~ q CR primary spectrum. These observations are encouraging in that they unambiguously 
confirm the breaks. However, they rule out traditional DSA models based on a single power law with an exponential cutoff. 

n. MECHANISM FOR THE SPECTRAL BREAK 

The physics of the spectral break considered here is very simple. When a SNR shock approaches a molecular cloud (MC) 
or a pre-supernova swept-up shell, confinement of accelerated particles deteriorates. Due to the particle interaction with mag- 
netic turbulence, confinement generally requires scales similar to the particle gyroradius iflll IT2I1 . However, strong ion-neutral 
collisions substantially enhance the role of particle pitch angle in wave-particle interaction. While the waves are in a strongly 
ionized (closer to the shock) medium they propagate freely in a broad frequency range at the Alfven speed Va = B/^/Anpi with 
the frequencies co = Wa- Here k is the wave number (assumed parallel to the local field B) and p,- is the ion mass density. As 
long as the Alfven wave frequency is higher than the ion-neutral collision frequency v,„, the waves are weakly damped. When, 
on the other hand, the ion-neutral collision frequency is higher (deeper into the cloud) , neutrals are entrained by the oscillating 
plasma and the Alfven waves are also able to propagate, albeit with a factor ypT/Po < 1 lower speed, where po is the neutral 
density. The propagation speed reduction occurs because every ion is now "loaded" with po/p; neutrals. In between these two 
regimes Alfven waves are heavily damped and even disappear altogether for sufficiently small p,/po -C 0.1. The evanescence 
wave number range is then bounded by k\ = v,„/2Va and kj = 2yJ pi/ PoVm/VA- These phenomena have been studied in detail 
in I33II34I1 . and specifically in the context of the DSA in J35l - [37ll . Now we turn to their impact on the particle confinement and 
emissivity. 

In the frame work of a quasilinear wave-particle interaction the wave number k is approximately related to the parallel (to 
the magnetic field) component of the particle momentum by the cyclotron resonance condition kp^/m = ±co c , where the 
(non-relativistic) gyro-frequency co c = eB/mc. Note that the appearance of = pji, where fx is the cosine of the pitch angle 
(see FigQ]), instead of the often used "sharpened" ll38tl resonance condition kp/m = ±co c is absolutely critical for the break 
mechanism. 

The frequency range where the waves cannot propagate may be conveniently translated into the parallel momentum range 

pi < \p\\\ < p 2 , (1) 

with 

pi = 2VAtn(i) c /Vi n , P%=-^ Vpo/Pi > Pi- (2) 

That a spectral break must form at the photon energy corresponding to the particle momentum p =p\ = pbr, can be readily 
understood from FigQ] The 'dead zones' pi < \p^ | < p2 imply that particles with \p^ | > p\ do not turn around (while moving 
along the magnetic field) and escape from the region of CR-dense gas collisions at a p\\/p fraction of the speed of light. More 
specifically, particles with pi < \pn | < p2 escape because they are not scattered, whereas particles with |/?|| | > p%, while being 
scattered, maintain the sign of pn, as they cannot jump over the gap pi < \p« | < pi and so escape as well. An exception to this 
are particles with sufficiently large p± that can be mirrored across the gap or overcome it via the resonance broadening. We 
return to this possibility later. 

The break can also be explained in terms of the confinement times of different groups of particles introduced above. Particles 
with |/»|| J > pi spend only short time T esc ~ L c /c (where L c is the size of the clump) inside the gas clumps. They propagate 
ballistically and their scattering time is assumed to be infinite, as there are no waves they can interact with resonantly (pi < 
\p\\ \ < pi) or they cannot change their propagation direction ( \pn | > p2). Particles with \p« | < pi are, on the contrary, scattered 
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Figure 1: Momentum space of accelerated protons. Particle scattering zones on the yp\\,p±J- plane of momentum space. Protons in the 

stripes p\ < |p||| < p2 are not scattered by waves (see text). Therefore, particles from the domains |p||| > p2 maintain their propagation 
direction and promptly escape from the dense gas region. 

intensively in pitch angle, they frequently change their direction, and so sit in the clump for T con f ~ L% / 'k ~ L^. / 'c 1 ; T/ SC . Here T sc is 
their pitch-angle scattering time and K is the associated diffusion coefficient. Not only T con f ^> T esc is required, i.e., T sc <C L c /c, 
but also T con f > L c /U s h, which means that the shock precursor is shorter than the clump k / £/ sn ;$ Lqr < L c (here £/ sn is the shock 
velocity, and Lqr is the thickness of the CR front near the shock). The last condition ensures that particles with p« > p\ that 
escape through the clump after having entered it from the shock side, will not interact with the shock after they exit through the 
opposite side of the clump, thus escaping upstream, Figj2] The reason for that is a low level of Alfven wave turbulence ahead of 
the CR precursor. We also assume that the ambient magnetic field does not deviate strongly from the shock normal, in order to 
allow these particles to escape through the far side of the clump. 

While particles with p> pi escape from the regions of enhanced gamma radiation (high gas density), an initially isotropic dis- 
tribution of accelerated particles is maintained only in a slab in momentum space \p^ | < pi and becomes thus highly anisotropic 
(a 'pancake' distribution). What matters for the integral emission, however, is a locally isotropic component / of this new proton 
distribution. It can be introduced by re-averaging the 'pancake' (|/»|| | < p\) distribution in pitch angle, f(p) = f f (p,n)djX, 
and is readily obtained assuming that particles remaining in the dense gas (those with \p^ | < p\) maintain the flat pitch-angle 
distribution, i.e. 



7M- //»«* = i ( ^W W '^: (3) 

where fo (p) is the initial (isotropic) distribution function and fj,i = min{p\/p, 1}. Thus, the slope of the particle momentum 
distribution becomes steeper by exactly one power above p = p\ = p\, T . In particular, any power-law distribution °< p~ q , upon 
entering an MC, turns into p~ q ~ l at p > p^, and preserves its form at p < p\, T . 

Note that the broken power-law spectrum can only be maintained if the filling factor / gas of the dense gas with the significant 
wave evanescence interval (p\ ,p2) is relatively small, / gas <C 1, so that the overall particle confinement and thus the acceleration 
are not strongly affected. If, on the contrary, / gas ~ 1, the resonant particles would leak into the (pi,p2) gap and escape from the 
accelerator in large amounts, thus suppressing the acceleration. We discuss further limitations of the mechanism in SeclVIIbut 
we note here that particles with sufficiently high momenta p > P2B0 / 8B, where 8B/Bois the effective mirror ratio of magnetic 
perturbations, can "jump" over the gap . The primary p~ q slope should then be restored for such particles. Recent MAGIC 
observations of the SNR W51C 0261 1310 indeed point at such spectrum recovery at higher energies. It should also be noted, that 
the Aq — 1 break index is a limiting case of identical gas clumps. The integrated emission from an ensemble of clumps with 
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Figure 2: SNR shock propagating into dense gas environment. The filling factor of the gas clumps is small, while some of them may be 
larger than the thickness of the CR layer near the shock front. 



different p\ and p2 should result in a more complex spectrum. 

III. BREAK MOMENTUM 

While the one power spectral break in the pitch-angle averaged particle distribution seems to be a robust environmental 
signature of a weakly ionized medium into which the accelerated particles propagate, the break momentum remains uncertain. 
According to eq.©, pb r (= pi) depends on the magnetic field strength and ion density as well as on the frequency of ion-neutral 
collisions, v,„ = no (oV) ■ Here (<JV) is the product of the collision cross-section and collision velocity averaged over the thermal 
distribution. Using an approximation of ll36ll39ll for (ffV), Pbr can be estimated as 

p br /mc~ lOB^XV 172 - (4) 

Here the gas temperature Ta, is measured in the units of 10 4 /T, magnetic field -in microgauss, no and n, (number densities 
corresponding to the neutral/ion mass densities po and p,) -in cm~ 3 . Note that the numerical coefficient in the last expression 
may vary depending on the average ion and neutral masses and can be higher by a factor of a few ll33i l4fjll than the estimate 
in eq.© suggests. The remaining quantities in the last formula are also known too poorly to make an accurate independent 
prediction of the position of the break in the gamma ray emission region. Those are the regions near the blast wave where 
complicated physical processes unfold. They include particle acceleration, strong MHD turbulence (driven by particles and their 
interaction with ambient gas inhomogeneities), gas ionization by shock generated UV photons, turbulent plasma heating and 
even evaporation of magnetic cloudlets 11391 |4U 14211 . Also important may be the ionization by the low energy CRs accelerated 
at the blast wave. However, as their diffusion length is shorter than that of the particles with p > pt, n we may assume that 
they do not reach the MC. Pre-ionization by the UV photons can also be ignored for the column density N > 10 l9 cm~ 2 ahead 
of the shock beyond which they are absorbed ll29ll . The authors or the Ref. [29], using the earlier data from ll43l have also 
analyzed the parameters involved in eq.(|4j and found the above estimate of pi, r to be in a good agreement with the spectral 
break position measured by the Fermi LAT. Nevertheless, we may run the argument in reverse and use the Fermi observations 
|01 of the gamma-ray spectrum of SNR W44 to determine the break momentum in the parent particle spectrum and constrain 
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Figure 3: Spectra of accelerated protons and electrons. The both particle distributions are calculated for a weakly modified shock and are 
shown in momentum normalization (/ (p) is steeper by two powers than the spectra in energy normalization, used in the text). Both spectra are 
multiplied by p A , so that the test particle distribution is flat. Shock parameters: acoustic Mach number M = 30, shock velocity V s /c = 10 -3 , the 
break momentum pi, r ~ line. Shock pre-compression (flow compression across the CR precursor) R=1.8, injection parameter v ~ 0. 1 [defined 

as v = (4^/3) {tru? lpVj\ (pinj/mc) f(pinj), with p and V s being the ambient gas density and the shock speed, respectively]; injection 
momentum p m jjmc ~ 1.4- 10~ 3 . 

the parameters in eq.®. Since we also know the amount of the slope variation Aq, we can calculate the full spectrum up to the 
cut-off energy. 

It should also be noted that in reality the break at p = p^ is not infinitely sharp for the following reasons. The break momentum 
may change in space due to variations of the gas parameters (eq.|4]), the resonance broadening [44, 45] near p = pi = Pbv ( so 
that particles with /? > pi are still scattered, albeit weakly) and other factors, such as the contribution of small gas clumps with 
L c <C Lcr, Figj2] The small clumps are submerged in the CR front and the CRs that escape from them are readily replenished. 
Note that this effect may decrease the break index Aq. However, the conversion of the parent proton spectrum into the observable 
gamma emission introduces a significant smoothing of the break, so that even a sharply broken proton spectrum produces a 
smooth gamma spectrum. It provides an excellent fit to the Fermi gamma data without an ad hoc proton break smoothing 
adopted by the Fermi-team [4] to fit the data. This will be seen from our calculation of the gamma emission based on the sharp 
proton spectral break in SecJVIbelow. 
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IV. PARTICLE SPECTRA 

To calculate the particle spectra, we need to determine the degree of nonlinear modification of the shock structure. In principle, 
it can be calculated consistently, given the shock parameters and the particle maximum momentum, p max - In the case of a broken 
spectrum, pt, r likely plays the role of p„ mx , as a momentum where the dominant contribution to the pressure of accelerated 
particles comes from, thus setting the scale of the modified shock precursor. Note that in the conventional nonlinear (NL) 
acceleration theory, the cut-off momentum p max plays this role, because the nonlinear spectra are sufficiently flat so as to make 
the pressure diverge with momentum, unlike broken spectra. 

The break in the photon spectrum is observed at about 2 GeV, which places the break in the proton distribution at about 
pt, r ~ IGeV/c [4]. For the strength of the break Aq=l, the spectrum above it is clearly pressure converging, and we perform 
the calculation of the shock structure and the spectrum using this break momentum as the point of the maximum in the CR 
partial pressure. Note that outside of gas clumps the CR pressure may still come from higher momenta and the complete 
nonlinear calculation of the spectrum would require the filling factor of the gas clumps. However, once the break momentum is 
set, we can use an analytic approach 11131 14611 for a stationary nonlinear acceleration problem using pt, r as an input parameter. 

Apart from pj, r , the nonlinear solution depends on a number of other parameters, such as the injection rate of thermal particles 
into acceleration, Mach number, the precursor heating rate and the shock velocity V s . Of these parameters the latter is known 
reasonably well, V s ~ 300km/ s, the injection rate can be either calculated analytically for the parallel shock geometry 114711481 . 
or inferred from the simulations N49L while the other parameters are still difficult to ascertain. Fortunately, in sufficiently strong 
shocks the solution either stays close to the test particle (TP) solution (leaving the shock structure only weakly modified) or 
else it transitions to a strongly modified NL-solution regime. The TP regime typically establishes in cases of moderate Mach 
numbers, low injection rates and low p mclx (now probably closer to pt, r ), while the NL regime is unavoidable in the opposite part 
of the parameter space. 

In the TP regime the spectrum is close to a power-law with the spectral index 2 throughout the supra-thermal energy range. 
In the NL regime, however, the spectrum develops a concave form, starting from a softer spectrum at the injection energy, with 
the index q ~ [r s + 2)/(r s — 1) > 2, where r s < 4 is the sub-shock compression ratio. Then it hardens, primarily in the region 
p ~ mc, where both the partial pressure and diffusivity of protons change their momentum dependence. The slope reaches its 
minimum at the cut-off (break) energy, which, depending on the degree of nonlinearity, can be as low as 1 .5 or even somewhat 
lower if the cut-off is abrupt. The question now is into which of these two categories the W44 spectrum falls? Generally, 
in cases of low maximum (or, equivalently, low spectral break pt r < 10) momentum, the shock modification is weak, so the 
spectrum is more likely to be in a slightly nonlinear, almost TP regime. On the other hand, there is a putative indication from the 
electron radio emission that their spectrum may be close to q e s» 1.75, which could be the signature of a moderately nonlinear 
acceleration process. It should be remembered, however, that this is a global index across the W44 remnant. There are resolved 
bright filaments where a canonical a = —0.5 spectrum, corresponding precisely to the TP parent electron spectrum with q e = 2 
is observed j50Tl . Moreover, there are regions with the positive indices a < 0.4 which cannot be indicative of a DSA process 
without corrections for subsequent spectral transformations such as an absorption by thermal electrons l50il. These regions 
can very well contribute to the overall spectral hardening discussed above, mimicking the acceleration nonlinearity. Finally, 
secondary electrons give rise to the flattening of the radio spectrum as well I29I1 . 

If the accelerated protons and electrons respond to the turbulence similarly, which is almost certainly the case in the ultra- 
relativistic regime, their spectra should have similar slopes there (as long as the synchrotron losses are ignorable). In using 
the electron radio spectrum as a probe for the level of acceleration nonlinearity, the following two relations are useful. First, 
there is a relation between the electron energy and the radio frequency Vmhz. = 4.6 -B^E^y. The second, already mentioned 
relation, q e = 1 — 2a, links the spectral index of radio emission a (assuming the radio flux « v a ) and the spectral index of the 
parent electrons q e (assuming their energy spectrum °< E~ qe ). Once the global radio spectral index of W44, a ~ —0.37 iTsoh 
is generated by freshly accelerated electrons in the frequency range 74 < v < 10700 MHz, the electrons should maintain their 
modified spectrum over the energy range spanning more than one order of magnitude. For example, assuming Bp, ~ 70 [4], 
one sees that electrons must maintain an index q e 1.75 between 0.46 < E < 5.8 GeV. While the upper bound is acceptable 
given the spectral break proton energy inferred from the super GeV emission measured by the Fermi LAT, the lower end is 
rather uncomfortable, since the nonlinear hardening of both protons and electrons with the Bohm (or other similar for protons 
and electrons turbulent diffusivities) starts (slowly) only at the proton rest energy. The calculated nonlinear spectra are shown in 
Fig|3]for the both species. At and below lGeV, the electron spectrum stays close to the test particle solution, q e K,2, even though 
the proton spectrum may steepen there, as we mentioned above. The physical reason for this difference is that the electron mean 
free path falls off slower with decreasing momentum than that of the protons in the Mev-Gev momentum range so that electrons 
sample longer parts of the shock precursor with higher flow compression and thus develop a harder spectrum. 
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Figure 4: Gamma radiation spectra. Photon spectra resulting from n decay and calculated for two different parent proton spectra compared 
against the Fermi (circles) and Agile (squares) data. Solid line: a test particle acceleration regime with the spectral index q = 2 below the 
break and q = 3 above the break at p/, r = IGeVlc. Dashed line: a moderately nonlinear acceleration regime corresponding to the spectrum 
shown in Fig[3](g ~ 1.75 and q ~ 2.75 below and above the break, respectively). Cut-offs are placed at 300 GeV for TP- and 100 GeV, for 
NL-spectrum. Fermi and Agile data are adopted from [4, 9], respectively. Both curves are well within the error bars of Fermi and Agile (not 
shown for clarity), which, in turn, overlap [9]]. 



V. PHOTON SPECTRA 



The above considerations somewhat weaken the radio data as a probe for the slope of the electron and (more importantly) for 
the proton spectrum. Therefore, the exact degree of nonlinearity of the acceleration remains unknown and we consider both the 
TP and weakly NL regimes in our calculations of the photon spectra, generated in p — p collisions. Specifically, we calculate 
the 7T° production rate and the gamma-ray emissivity. In so doing, we adopt numerical recipe described in detail in IdU 1531 . 
The physical processes behind these calculations are (i) collisions of accelerated protons with the protons of the ambient gas 
resulting in the following spectrum of 7T°-mesons: 

F pp (E K ) = 4%N pg J da ^*f p) J P (E P ) dE p 

where, N pg is the number density of protons in the gas, do /dE K is the differential cross section for the 7T° production in collisions 
between accelerated protons of energy E p and gas protons, J p is the flux of accelerated protons, E K is the energy of n° mesons; 
(ii) decay of n° resulting in the gamma emission spectrum 
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J y/E*-nfy* 

E r +mic 4 /4E Y 

where m n is the pion rest mass. 

The results are shown in Figj4] The best fit to the Fermi and Agile data is provided by a TP energy distribution («; E~ 2 ) below 
Pb r — 7GeV/c with the spectrum steepening by exactly one power above it. The spectrum steepening is perfectly consistent 
with the proton partial escape described above (with no parameters involved) and shown in Fig[T| For comparison, a weakly NL 
spectrum, shown in Figj3] is also used for these calculations (dashed line in FigHJ, but its fit would require a somewhat stronger 
break (Aq > 1) or a low momentum cut-off, Figj4] i.e. at least one additional free parameter. We will discuss the options in the 
next section. 



VI. SUMMARY AND OUTLOOK 

To summarize the results, the mechanism for a break in the spectrum of shock accelerated protons suggested in lEM l30ll is 
in excellent agreement with the recent [4] Fermi LAT and Agile [9] observations of the SNR W44. The observed gamma ray 
spectrum most likely results from the decay of 7C°-mesons which are born in p — p collisions of shock accelerated protons with 
an ambient dense gas. The parent proton spectrum is best represented by a classical test particle power law oc E~ 2 , steepening to 
at Eb r s» IGeV due to deteriorated particle confinement caused by the ion-neutral collisions and the resultant Alfven wave 
evanescence. The position of the break momentum in the particle spectrum may be estimated using eq.©, or conversely, the 
combination of parameters involved in this estimate can be inferred from the measured break momentum. The cut-off momentum 
is not constrained in this scenario. 

An alternative explanation, based on a different mechanism of the break, associated with the change of the particle transport 
in the CR shock precursor 05311 is also possible but is less definitive in the spectrum slope variation Aq across the break (see also 
l29ll for the most recent alternative suggestions). In addition, the mechanism [53] would imply a considerable nonlinearity, i.e. 
a stronger CR shock precompression than that suggested by the radio observation of accelerated electrons and the inferred 100 
GeV proton upper cutoff (see below). Still alternatively, assuming the "environmental" break mechanism is at work, i.e. Aq=l, 
but the shock structure is somewhat modified, we arrive at the E~ 1J5 spectrum below the break (as the radio observations may 
suggest for the electrons), and E~ 2J5 above the break. A fit to the data is marginally possible, but it would require a relatively 
low cut-off momentum at about 100 GeV/c. This possibility may be supported or ruled out once the data (upper limit) around 
this energy become available. 

As we noted, particle escape from the MC can quench the acceleration process 05411 . This would certainly be the case if the MC 
were filling the entire shock precursor. However, MCs are known to be clumpy B55I458I1 . and fill only a small fraction (< 1-2%) 
of the precursor. In this case the acceleration process continues largely unimpeded (apart from the spectrum steepening) but the 
accelerated protons illuminate the 'cloudlets' and make them visible in y-rays due to the high density target material. Another 
concern is a faint or even lacking x-ray emission that seems to be inconsistent with shocks impacting dense surroundings. This 
issue has been recently dealt with in, e.g., ll58il59ll . Large clumps survive the shock passage as it stalls inside them and no strong 
heating occurs [|58|] . 

The most robust and attractive aspect of the suggested mechanism for the spectral break is the exact Aq = 1 variation of 
the spectral index. Indeed, this change in the spectral slope is due to the reduction of the number of degrees of freedom of 
particles caused by the resonant wave evanescence and it does not depend on any parameters. In a combination with the test 
particle regime operating below the break, which is physically suggested by the low values of the break and upper cut-off 
momenta, the mechanism provides a very good fit to the Fermi LAT and Agile data with no free parameters for the SNR W44 
and probably for W51C. From a number of physically different types of spectral breaks suggested [53, 60-62], namely the 
current, "environmental" mechanism appears to be plausible where a dense target gas is present which is also required for the 
efficient n° production. However, observations of some other remnants in the dense gas environments, such as W28 and IC443 
10, 0] indicate weaker breaks, Aq = 0.6 — 0.7 which may either require a different mechanism for the break or a narrower wave 
evanescence gap Ap = pi — p\ (higher ionization rate). The predominance of small clumps with L c <C Zcr m a MC will also 
reduce Aq. 

Generally, spectral breaks offer a natural resolution to the well known but puzzling trend of the nonlinear (i.e. supposedly 
improved) DSA theory to develop spectra which are considerably harder than a simple test particle spectrum, thus becoming 
even less consistent with the bulk of observations 063116411 . However, the nonlinear spectrum - i.e., diverging in energy- exhausts 
the shock energy available for the acceleration as the cut-off momentum grows, so that a broken spectrum should form N53ll6lll . 
Broken spectra are now commonly observed and the old paradigm of a single power-law with an exponential upper cut-off 
is maladapted to the recent, revolutionarily improved observations [0,0]. Note, that the spectrum of the RX J17 13.7-3946 
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^ is also consistent with the environmental break mechanism presumably operating in W44 surrounding but with a higher 
Pbr ~ 10 3 GeV jc and thus with stronger acceleration nonlineari ty p 8ll. However, it is difficult to make the case for hadronic 
origin of the gamma-ray emission of the RX J17 13.7-3946 J3l|65l|6a]. The fundamental role of the W44 remnant for the problem 
of CR origin is that this particular remnant seems to rule out contaminating electron emission due to Bremsstrahlung and inverse 
Compton scattering (j, |29J thus favoring the hadronic origin of the gamma emission and bolstering the case for the SNR origin 
of galactic CRs 
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